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Abstract 

The paper considers a linear matrix inequality (LMI) that depends on a parameter varying in a 
compact topological space. It turns out that if a strict LMI continuously depends on a parameter 
and is feasible for any value of that parameter, then it has a solution which continuously depends 
on the parameter. The result holds true for LMIs that arise in S-procedure and Yakubovich 
lemma. It is shown that the LMI which is polynomially dependent on a vector of parameters 
can be reduced to a parameter-independent LMI of a higher dimension. The result is based 
on the recent generalization of Yakubovich lemma proposed by Iwasaki and Hara and another 
generalization formulated in this paper. The problem of positivity verification for a non-SOS 
polynomial of two variables is considered as an example. To illustrate control applications, a 
method of parameter-dependent Lyapunov function construction is proposed for nonlinear sys- 
tems with parametric uncertainty. 
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1. Introduction 

The S-procedure losslessness theorem and the Kalman-Yakubovich-Popov (KYP) lemma are 
important mathematical tools of modern control theory. Both statements deal with inequal- 
ities of certain quadratic forms. We consider a situation when these quadratic forms depend 
continuously on a parameter that varies in a compact topological space. 

First, we study solutions of a parameter-dependent linear matrix inequality (LMI). It turns 
out that if a strict LMI continuously depends on a parameter and is feasible, then it has a 
solution that continuously depends on the parameter. From this general statement it follows 
that if S-procedure for strict inequalities is lossless for any value of the parameter, then the 
Lagrange multipliers can be chosen as continuous functions of the parameter. 

We also consider parameter-dependent generalizations of one statement of KYP lemma, 
namely the equivalence of strict frequency- domain inequality and the strict linear matrix in- 
equality. (This result was first proven by V.A. Yakubovich [1] and was called Yakubovich 
lemma by R. Kalman [2] and S. Lefschetz [3].) 

It is proven that if the matrices in the frequency- domain inequality continuously depend 
on the parameter, then there exists a solution of the LMI, which continuously depends on the 
parameter. The result also holds true for the generalized KYP lemma proposed by T. Iwasaki 
and S. Hara [I] as well as for a new version of KYP lemma with matrix frequency- domain 
inequality. The different result concerning a parameter-dependent version of KYP lemma was 
obtained by A. L. Likhtarnikov [5]. In [5] solutions of Lur'e equation are considered that are 
analytic functions of a parameter. 

Using parameter-dependent Yakubovich lemma we show that if an LMI polynomially depends 
on the parameters, then it can be transformed into an LMI of a larger dimension that does not 
depend on these parameters. As an illustration of this result, we consider minimization of a 
polynomial of several variables in a bounded domain. The proposed method can be used to 
verify the positivity of polynomials that cannot be represented as a sum of squares. 

The results are also applied to the construction of a parameter-dependent Lyapunov function 
for stability analysis of nonlinear systems with parametric uncertainty. 

2. Parameter-dependent LMI 

Consider the parameter-dependent LMI in h G R 

L(p,h)>0 (1) 

where p G T> is a parameter, L : T> x R' — > HM n is affine with respect to h. Hereafter HM„ 
denotes the space of Hermitian matrices of dimension n. 

Theorem 1 Suppose that V is a compact topological space, the map L is continuous and is affine 
with respect to h, and that LMI (TJP is feasible for any p G T>\ then there exists a continuous 
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function h(.) : T> — > M. 1 that satisfies inequality 

L(p,h(p)) > Vp G V. (2) 
The set of all continuous h(.) satisfying is open in C(Z>,R ). 

Hereafter C(X, y) denotes the normed space of continuous functions from the topological space 
X to the normed vector space y . 

Corollary 1 If T is a dense subset of C(V,M. 1 ), then there exists h(.) G T satisfying (TJ|). 

Corollary 2 1/Dc R , then there exists a polynomial h(p) satisfying (d)). 

Suppose that h{p) is a polynomial of a vector variable p G R fc . Denote the vector of coefficients 
of h(p) by h G R fc . Let L : V x R k -> HM n be a map that is defined by substitution of the 
polynomial h(p) into L. From corollary [2] it follows that 

L(p, h)>0 Vp G V. (3) 

This way the construction of the parameter-dependent solution h{p) of ([1]) is reduced to a 
search of a constant vector /i that satisfies the system of inequalities ([3]) indexed by p. The 
most interesting case is the one when this system includes infinite number of LMIs. In section 7 
we show that under certain assumptions on T> this system can be reduced to a single finite- 
dimensional LMI. 

3. Parameter-dependent S-procedure 

Consider two maps i*o : T> x X — > R and F : T> x X —>■ M 1 , where T> is a compact space, 
^ is either C n or R n . Let /C be closed convex proper cone in R'. The dual cone is denoted by 
/C + = {z G R z |(2, y) > Wy G /C}. Consider two statements: 

I. Vp G X> F(p, x) G /C =^ F (p, x) > 0. 

II. 3h:V^JC + such that 

F (p,a;) > (h(p),F(p,x)} VpeV lX e X. (4) 

It is clear that II implies I. In the case when parameter p is absent, the operation of re- 
placement of the condition I by the condition II is known as S-procedure. In the presence of 
parameter we call this method a parameter-dependent S-procedure. The S-procedure is said to 
be lossless if I is equivalent to II. 

Theorem 2 Suppose that Fq,F are continuous and are quadratic forms of x. If parameter- 
dependent S-procedure is lossless then there exists a continuous function h(.) that satisfies (TJJj. 
The set of all continuous h(.) satisfying (fj]j is open in C(T>,W). 
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4. Generalized Yakubovich Lemma 



Let us begin from the formulation of the classical result. 

Theorem 3 (Yakubovich 1962, 1964) Let A G O xn * , B G (O xnu , G G MM nx+riu . The 
following statements are equivalent: 

1. The inequality 

':)* g (:)>° < 8 ' 

is fulfilled for all x G C nxXl , u G C nuXl , |x| + \u\ ^ 0, u G [— oo, +oo] snc/i iwx = Ax + Bu. 

2. There exists H G HM ni that satisfies the LMI 

<«> 

Hereafter i is the imaginary unit. 

This theorem was first proven by V.A. Yakubovich in [I] for the single-input system (i.e. for 
n u = 1) and in [6] for the multi-input system. Our formulation follows |7J. Now the statement 
of the Theorem [3] is known as a part of Kalman- Yakubovich- Popov lemma [TJ [2j [6j [8] . In our 
opinion, it is more convenient for reference purposes, and is more correct from the historical 
point of view, to refer to Theorem [3] as Yakubovich Lemma. We call this statement a lemma 
to point out its connection with Kalman- Yakubovich-Popov Lemma. Additional information 
about history of Kalman- Yakubovich-Popov Lemma can be found in [9]. 

The theorem below is proven in [1]. It is a generalization of Yakubovich Lemma. The paper 
[I] is based on generalizations of Yakubovich Lemma that were proposed in [TU1 [TT] . 

Let 6 = 6* = ( 6 Q n ° Q 12 ) , det 6 < 0. Define the curve r e = {A G C | (A, 1)6(A, 1)* = 0}, 

V #21 t/ 22 / 

and the domain Q@ = {A G C | (A, 1)0 (A, 1)* > 0}. Hereafter C is the closed complex plane 
C = CU {oo}. Choosing different matrices B it is possible to define any circle or straight line on 

C. To illustrate, let us present some examples: ® = ^j* g^>re = iMU {oo},fie = {ReA > 

0}U{oc}; 6= ( •? " o 1 ^,r = MU{oo},fi = {ImA>O}U{oo};0= ( ~ 1 ° J , T e = 

{|A|=r},fi = {|A|<r}. 

Let M,N G C nxn % S G iM„ z , H G IM„. Define the generalized Lyapunov opera- 
tor A m ,n,b(S) = (M,N)(e <g> S)(M, N)*. The adjoint operator takes the form A' MN e (H) = 
(M*, N*)(Q T ®H)(M*, N*)*. Hereafter M 1 ®M 2 is a Kronecker product of matrices M x and M 2 . 

Theorem 4 (Iwasaki and Hara, 2005) Let M, N G C nxriz , G G HM„ z . Suppose the inter- 
section Tq 1 fl Qq 2 includes more than one point and n z > n; then the following statements are 
equivalent: 

1. The inequality 

z*Gz > (7) 
is fulfilled for all z G O xl , z ^ 0, A G r 6l n fi 02 such that (AY - M)z = 0. 



4 



2. There exist Hx,H 2 E HM n , H 2 > such that 



(8) 



5. Yakubovich Lemma for matrix frequency-domain inequality 

Let M,N,Qi,i = 1,2 be defined as in previous section. Consider matrices Gy G C" zXnz , 
Gij = G* j{ , i,j = 1,. . . ,m. 

Theorem 5 Suppose the intersection Tq 1 n fie 2 includes more than one point, n z — n+ 1, and 
rank(AiV — M) = n /or a// A G 1^ n ^e 2 ; ^en i/ie following statements are equivalent: 
1. T/ie matrix inequality 



*G\ m z \ 



> 



(9) 



y 2? G m x2^ ... z G mm z J 



is fulfilled for all z G O xl , z ^ 0, A G i\ n fi 2 suc/i ffcat (AiV - M> = 0. 
2. T/iere exzsi #i,# 2 G HM mn , H 2 > snc/j i/iai 



GmA ... G r 



(10) 



/ 



6. Parameter-dependent Yakubovich Lemma 



Let £> be a topological space. Consider maps Gij(.) G C(X>, C nzXriz ), (Gy(p) = G 
£>, i,j = l,...,m), M(.),iV(.) G C(£>,C nxnz ), 9;(.) G C(P,BM 2 ), (det9;(p) <0Vp6P, i = 

/ G u (p) ... G lm (p) \ 
1,2). Denote G(p) = ; ■-. ; 

\ G TO i(p) . . . G mm {p) j 

Theorem 6 Suppose V is a compact space, the intersection ^&i(p) D ^e 2 (p) includes more than 
one point for all p G P, and either m = 1, n z > n orm > l,n z = n+1, rank(AiV(p) — M(p)) = n 
/or a// A G Te^p) PI fie 2 (p) and a// p G T>; then the following statements are equivalent: 
1. Tne matrix inequality 



I **Gh(p); 



'Glm(p)« \ 



> 



(11) 



\ z*G ml (p)z ... z*G mm (p)z J 



is fulfilled for all p G £>, z G O xl , z ^ 0, A G r 6l (p) H ft 02 (p) suc/i too* (AiV(p) - M(p))z = 0. 
2. T/iere exzsi H 1 (.),H 2 {.) G C(D,iM mn ) stzc/j tfiat /or all p e V 

H 2 (p) > 0, G(p) > A' ImmHp ),i m m(p),e 1 (p)( H i(p)) + A L®M(p),/ m ®7v(p),e 2 (p)(#2(p))- (12) 
The set of continuous pairs (Hi(.),H 2 (.)) satisfying ( THij) zs open m C(P, HM mn x HM mn ). 
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7. Parameter-independent solutions of parameter-dependent LMI 



Consider the parameter-dependent LMI (JTJ. In contrast to section 2 now we are looking 
for a constant vector h that satisfies (CQ) for all p G T>. In section 2 it was shown that when 
Pel' the search of parameter-dependent solution of an LMI can be reduced to the search of 
constant one for a set of parameter-dependent LMIs with growing dimension of solution. Now 
we would like to show that in some cases the search of the parameter-independent solution of 
a parameter-dependent LMI can be reduced to the solution of a set of parameter-independent 
LMIs with growing dimension. 

Suppose that L is polynomial of p and the set V is given by 



V = {p g M lx >i < Pl < h, ai (p 1 - 1 ) <pi< b^p 1 - 1 ), i = 2,...,k}, 



(13) 



where p = (pi, . . . ,Pk), P % = (pi, ■ ■ ■ ,Pi) G R lxj , a± < b\, ft, b iy i — 2, . . . , k, are polynomials and 
< hip 1 - 1 ) for all p G V. 
Let us describe the procedure for construction of a parameter-independent LMI, each solution 
of which defines the vector h satisfying ([I]) for all p G T>. 

de gpk L +1 j j ((p k ) = (pf, . . . ,p k ,l) G R lx ( d+1 ); then L(p, h) can be represented as 



Define d — |_- 



L(p,h) 



( (Gu^-WC ... {G lm (p k -\h)(* \ 

yCGmlip^^h^* ... CGmm{p k -\h)C J 



where the matrices Gij(p ,h), i,j = l,...,m, are polynomials of p and are affine with 
respect to h. For all i,j — l,...,m, G i j(p k ' 1 , h) = Gj i (p k ^ 1 , h)*. 
Define matrices M = (I d , 0), N = (0, I d ) G R dx ( d+1 \ 



-i 

1 



,e 2 (p 



1 



\M P k - l ) + b k {p k - 1 )) 



\{a k {p k 1 )+b k (p k ~ 1 )) -a^p^bkip 



k-l\ 



Then 



{p k N - M)z = 0^3ceC: z = c(p d k , ...,p k , 1) 

r 01 n fie 3 (p*-i) = kb* -1 ),^ -1 )]- 

Thus (P) is fulfilled for all p G £> iff 



> 



(14) 



\ z*G ml (p k \h)z ... z*G mm (p k \h)z J 



for all p*- 1 G P*- 1 , 2 G C( d+1 > xl , z ^ 0,p fc G r 6l n ^e 2 (p*-i), satisfying (p fc iV - iW> = 0. Here 
V k - X = {q G R^^-^lai < ft < 6i, a^- 1 ) <qi< i = 2, . . . , k - 1}. 

Let us introduce matrix polynomials i/i(p fe_1 ), H 2 (p k ~ 1 ) G C('D fe_1 , HM^) . Consider in- 
equalities 



G(p k -\h)>A' 



e^Hiip*- 1 )) + A', 



-m ®M,I m <S>N,Q 2 (q) 



(H^- 1 )), H 2 {p^) > 0. (15) 
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where 

/ G u (p k -\h) ... G lm (p k ~\h) \ 

G(p k ~\h) = : ■.. : EUM m(d+1) . 

\ G ml (p k ~\h) ... G mm (p k -\h) J 

From Theorem 6 it follows that if the vector h and the pair of polynomials Hi(p k ~ l ), H 2 (p k ~ 1 ) 
satisfies (fT5j) for all p fc_1 G £> fc_1 , then h satisfies ([Tj) for all p G V. 

Let haux be the vector of coefficients of polynomials Hi(p k ~ 1 ), H 2 (jp k ~ l )\ then inequalities 
( TT5|) can be rewritten as follows: 

Ri, i (jf-\h,h£> x )>0, i = i,2, (16) 

where Ri,i{p k ~\ h, h£L) = G{p k -\ h) - A/ m ®M,/ m ®Jv > G 1 (^i(P fe-1 )) - ^^^^^^^^(^(p^ 1 )), 
R^-^Kh^l) = H 2 (p k - V ). The maps it^j, z = 1,2, are polynomials of ^ 1 and are affine 
with respect to the joint vector (h, haux)- 

Applying the same procedure to each inequality in (fT5"j) . we obtain the system of four in- 
equalities 

R 2 ,i(p k ~ 2 , h, h£l) > 0, i = 1, 2, 3, 4, (17) 

(2) 

where haux is the vector of parameters of all polynomials, that were introduced on two steps of 
procedure, the maps R 2 ,i, i = 1,2,3,4, are polynomials of p k ~ 2 and are affine with respect to 

(2) 

the joint vector (h,h a u X ). 

Repeating the procedure fc-times, we obtain the system of parameter-independent LMIs that 
has the form 

Rk,i(h, h aux ) > 0, z = 1, . . . , 2 fc , (18) 

where h aux is the vector of parameters of all introduced polynomials, the maps Rk,%, i = 1, ■ ■ ■ ,2 k , 
are affine with respect to the joint vector (h, h aU x)- 
For convenience we rewrite (|T8|) as a single LMI 

R{h,h aux )>0, (19) 

where R is block-diagonal matrix with blocks Rk,i, i = 1, . . . ,2 k , h aux G M. 1r . 

Let 1Z be the set of all affine maps R that can be constructed from the map L, using described 
recursive procedure. 

Theorem 7 Suppose that L is polynomial of p, D is given by iH~3\) ; then 

{heR l \ L{ P , h)>0V P eV}= J {h e R l | 3h aux e R lR : R(h, h aux ) > 0}. 

Rell 

8. Example. Minimization of a polynomial in a bounded domain 

Let g(p) be the polynomial of p G M lxfc . Consider the problem of evaluation of the polynomial 
g minimum in a domain T> given by (|13p . 
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Define L(p, h) = g(p) — h, h G M. It is clear that any h satisfying ([I]) for all p G T> is a lower 
bound for g ml). Moreover from Theorem 7 it follows that 



sup sup{/i | 3/w G R lR : /w) > 0} = mm{g(p) \ p G £>}. 



Ren 



(20) 



Taking into account (TSUI) , we can say that minimization of the polynomial g in a domain X> is 
reduced to a standard LMI optimization problem. 

To illustrate the method we apply this approach to a polynomial of two variables g(x, y). Let 
V = {(x, y) | a < x < b, c(x) < y < d(x)}, where a < b, c(x), d(x) are polynomials, c(x) < d(x) 
for all x G [a, b}. Consider the parameter-dependent LMI in h G 1R 



g(x, y) — h> 0, Vx, y EV 



(21) 



Let deg y g be the degree of the polynomial g with respect to y. Denote n y = |_ de ^ J, 
y = (y"*, . . . ,y,l). Then #(x, y) — h = YG(x, h)Y*, where 
/ Gi,i(cc) ... Gx >n+1 (x) \ 



G(x, h) 



G HM^+i, Gij(x),i,j = 1, . . . ,n y + 1, are 



\ C^+^^x) ... G n;/+ i inj/+ i(a;) — /i y 
polynomials of x. 

Define matrices M y = (I ny ,0), N y = (0, I Uy ) G M n « x K +1 ). Then 



(yN y -M y )Z = 



(22) 



implies Z = cY T , c G C. 
Define matrices 



6i 



-i 

1 



,e 2 (x) 



-1 ±{c{x) + d{x)) 

;(c(x) + d(x)) —c(x)d(x) 



Then Tq-^ fl fle 2 (x) = [c(x), d(x)]. Thus (1211) is equivalent to Z*G(y)Z > Vx G [a, 6], y G 
To! H fie 2 (x) and satisfying f l22l . 

By Theorem 6 it follows that the latter condition is fulfilled iff there exist polynomial matrices 
H G1 (x), H G2 (x) G C(f,iM % ) satisfying the LMIs 



H G2 (x) > 0, G(x,h) > K' MytNyjQi {H G1 (x)) + M My ,N y ,@ 2 {x) (H G 2(X)) 



Let n be the maximal degree of x in (T23]) . Put n x = ; X 



lx" 



(23) 

x, 1). Let us 



denote the vectors of the coefficients of the polynomials H G1 (x), H G2 (x) by H G i, H G2 respectively 
and define matrices E(H G2 ,x) = H G2 (x), F(H G1 ,H G2 ,x,h) = G(x,h) - A' MyjNvjei (H G1 (x)) + 
K' My N ^ e ^ x JH G2 (x)) . Then LMIs ff23|) can be rewritten as follows 



/ XE hl (H G2 )X* 



E(H G2 ,x) 



XE 1>ny (H G2 )X* \ 



> 0, 



(24) 



\ XE nyjl (H G2 )X* ... XE ny ^ ny (H G2 )X* J 
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/ XF lt i(H G1 ,H G2 )X* 



F(H G1 ,H G2 ,x,h) 



XF ln+1 (H G1 ,H G2 )X* \ 



> 0. 
(25) 



\ XF ny+ i )1 (HGi,H G 2)X* ... XF ny+1>ny+1 (H G1 ,H G2 ,h)X* ) 
Define matrices M x = (i^O), N x = (0, J n J G R««x(«»+i). Then 

(xiV x - M X )Z = (26) 
-1 5(0 + 6) 



implies Z = cX , c G C. Let 63 



2(0 + 6) — ab 



According to Theorem 5, < HM and ( |25i) are fulfilled iff there exist Hei,He2 G HM ni „ s , 
fl>i, ifp2 e HM n:c ( ns+1 ) that satisfy the LMIs 



H E2 > 0, H F2 > 0, 



/ F U (H G1 , H G2I 



G2j 



\ E nyA (H G2 ) ... E nyinv {H G2 ) J 

Fln y +l(H G i, H G2 ) \ 



A'j 



y ®M x ,i ny ®N x ,e 3 



, Qi (Hei)+ 

qSh E2 ), 



(27) 



V En y 



+n{H G i, H G2 



F ny +ln y +l(H G i, H G2 ,h) ) 



> K ny+1 ®M x ,I ny+1 <g>N x ,e 1 (.FfFl) + 



In this way (121]) is fulfilled iff there exists n x such that LMIs fl27|) are feasible. Evidently, 
system f j27|) can be written as a single LMI f|T9l) . Maximizing h over solutions of fl27j) and 
increasing n x we can estimate minimum of g in D with any desired accuracy. 

The proposed approach can be applied to verification of positivity of polynomials. As an 
illustration let us consider a numerical example. 

It is known that g{x,y) = x A y 2 + x 2 y 4 — 3x 2 y 2 + l + e > for all e > and all x,y. Besides for 
any e the polynomial g(x, y) cannot be represented as sum of squares (SOS) of polynomials [12J. 
It means that the well known SOS representation technique cannot be used directly to verify 
the positivity of g. 

It is easy to see that g(x, y) > 0, when \x\ > 2 or \y\ > 2. So, to verify the positivity of g it is 
sufficient to find positive lower bound of g(x,y) in V = {(x,y) \ \x\ < 2, \y\ < 2}. Calculations 
using Matlab® LMI toolbox showed that for n x = 2 the system of inequalities fl27|) has a solution 
with h > if e = 1CT 11 . This proves that x 4 y 2 + x 2 y 4 — 3x 2 y 2 + 1 + 1CT 11 > for all x, y. 

9. Construction of parameter-dependent Lyapunov function 

Consider the parameter-dependent nonlinear system 

x = A(p)x + B(p)u, x(0) = x , (28) 

u = ip(p,t,x), (29) 
where t > 0, x,x G « nxl ,u G M mxl , p is a parameter, p G V C M lxfc , A(.) : V -> « nxn , 



B(.):V 



\ ip : V x [0, +00) x 



P71X 1 



pmx 1 
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The nonlinearity ip satisfies the quadratic constraint 



(x*,u*)G(p)(x*,u*Y > (30) 

which is fulfilled for all x G W lXl , t > 0, p G V, u = tp(p, t, x). Here G(p) G HM m+n . 

Suppose that A, B, and G are polynomials of p, and V is defined by ( TIB"]) . Consider 
the parameter-dependent Lyapunov function candidate V(p, x) = x*H(p)x, where H{.) G 
C(V,MM n ). We are looking for H(.) that satisfies 

H(p) > Vp G V, (31) 

and 

^y(p,x(t)) <0 Vt>0, (32) 

for all solutions of (128]) satisfying (130|) . If fulfilled these conditions guarantee the asymptotic 
stability of the closed- loop system (J28l) . (!29l) for any such that (J30l) holds. 

Define matrices M (p) = (A(p), B(p)), N = (I n , 0) G R™ x ("+ m ) ; 6 = J V Then 

_ / H(p)A(p) + A*(p)H(p) H(p)B(p) \ 
A M(piNM n W)- { B*(p)H(p) J 

and condition fl32|) takes the form: 

Vp G P, V|x| + \u\ ^ (J30D implies (x*, ^)A' M(p)jiVie ( J ff(p))(x*, m*)* < (33) 

Using parameter-dependent S-procedure and Theorem 1 we can see that fl33l) is fulfilled iff there 
are a polynomial matrix H(p) and a polynomial rj(p) that satisfy 

A^^oC^Cp)) + V(P)G(P) > 0, n{p) > Vp G P. (34) 

Let /i G R ! be the vector of coefficients of polynomials H(p),t]{p)- Then (j34D can be written 
as a single parameter-dependent LMI jT]). Using the procedure described in Section 7 we can 
define the set of affine maps TZ. From Theorem 7 it follows that H(p) and nip) satisfying ( TSTI) 
and (1341) exists iff there is an affine map R G 7£ such that LMI (fT9l) is feasible. Any solution h 
of obtained in this way LMI (fT9l) defines the polynomial matrix H(p) that satisfies (!3T|) and (|32|) . 
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